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Abstract 

In this paper, we investigate the downlink resource 
allocation problem for on-demand CM services in HSR 
OFDMA systems with a cellular/info station integrated 
network architecture. The resource allocation problem is 
formulated by a two-stage optimization programming, 
which aims at maximizing the total reward of delivered 
services then minimizing the weight total number of 
cumulative discontinuity packets over the trip of the train. 
An equivalent one-stage programming is proposed to 
resolve the difficulty of multi-stage optimization. The result 
mixed integer programming (MIP) is NP-hard in general, 
thus we reformulate it as a sparse l 0 -minimization problem 
and then relax it to a linear programming (LP). Guided by 
the proposed optimization approach, we develop an efficient 
online resource allocation algorithm for practical systems, 
where a-priori knowledge of future service arrivals is not 
available. Finally, simulation results are provided to validate 
the proposed algorithms. 
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Introduction 

In recent years, the high-speed railway (HSR) has been 
developed rapidly as a fast, convenient and green 
public transportation system (J. Wang etc., 2014). 
There is a strong demand for various multimedia 
services such as online game and videophone, 
especially during a long time travelling (L. Tian etc., 
2012; S. F. Xu etc., 2014). However, the data 
transmission rate degrades significantly when the 
train moves with an extremely high speed, e.g. 360 
km/h, due to the serious Doppler shift, high 
penetration loss of signal into the carriage and varying 
channel conditions (Y. Zhao etc., 2012). Specifically, 
the passengers in the train will not enjoy the 
continuous multimedia (CM) services smoothly, since 


the CM services, such as IPTV and VoIP (O. B. Karimi 
etc., 2012), is transmitted as a continuous packet flow. 
Hence, how to deliver CM streaming effectively over 
limited resource and provide passengers with a high- 
quality end-user experience is an interesting and 
challenging problem in the HSR environment. 

To address these challenges, we consider a 
cellular/infostation integrated HSR network 
architecture (H. Liang and W. Zhuang, 2012) in this 
work. In the integrated network, a cellular network 
with seamless coverage is supposed to support control 
channels for service requests and acknowledgements, 
while infostations deployed along the rail line can 
provide high-rate, high-quality and efficient data 
transmission. Based on this architecture, an optimal 
resource allocation problem for on-demand data 
delivery is investigated (H. Liang and W. Zhuang, 
2012), which takes the intermittent network 
connectivity and multi-service demands into 
consideration. In (T. Chen etc., 2013) the total weight 
of received packets insteading of the fully completed 
requests was maximized. However, these works 
focused only in the data transmission integrity, with 
no attention to the continuity of data delivery, which 
will finally affect the playback continuity at passenger 
devices. For the best of our knowledge, the resource 
allocation problem that can jointly consider the 
integrity and continuity of on-demand data delivery in 
such an intermittently connected network has not been 
studied in the previous works. 

In this paper, we investigate the downlink resource 
allocation problem in HSR downlink OFDMA 
(Orthogonal Frequency Division Multiple Access) 
system with a cellular/infostation integrated network 
architecture. The inter-carrier interference (I Cl) due to 
the high moving speed is considered in characterizing 
the link capacity. Taking the train trajectory and CM 
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service demands into account, the resource allocation 
problem is formulated by a two-stage optimization 
programming, wherein the integrity of on-demand 
data delivery is considered in the first stage, while the 
continuity is considered in the second stage. Then we 
reformulate the two-stage optimization as a single- 
stage i Q -minimization problem, and it is further 
approximated by a linear programming (LP). 
Moreover, we propose an efficient online resource 
allocation policy based on the insights gained from the 
offline analysis. 

The rest of this paper is organized as follows. The 
system model and assumptions are described in 
Section 2, and the problem formulation is provided in 
Section 3. Section 4 presents the offline resource 
allocation algorithm and Section 5 proposes the online 
algorithm. Simulation results and discussions are 
given in Section 6, followed by the conclusions in 
Section 7. 

Notation: [x] + = max{x, 0}and [xj = max{ n G Z|n < x}. 
Z + represents the set of nonnegative integers. For 
given vector x, || x || 0 and || x \\ 1 denote the l 0 and 
norm, respectively. 

System Model 

As shown in Fig. 1, we consider a HSR downlink 
OFDMA system with a cellular /infos tati on integrated 
network architecture (H. Liang and W. Zhuang, 2012). 
H infestations with small coverage areas are deployed 
along the rail track to provide efficient data 
transmission, while the cellular network with seamless 
coverage is used for supporting the control channels 
over the region. A central controller (CC) is deployed 
to allocate the network radio resources for the whole 
network. The infestations can communicate with the 
vehicle station (VS) installed on the top of the train. 
Meanwhile, the VS is further connected to the access 
points (APs) inside the train. Assuming that the data 
transmission rate from the VS to passenger devices is 
sufficiently high, hence the data packet can be 
received successfully if it has been delivered to the VS. 
When the train is moving, the passengers send service 
requests from the VS to the content server (CS) 
through the cellular network. The requested data 
packets are then delivered from the CS to the VS via 
the infestations, and the VS will eventually forward 
the data to the passenger devices. For simplicity, we 
assume that the buffer space of VS is unlimited and 
the data service rate is equal to the end-user playback 
rate, which is assumed to be constant. 



Note that each high-speed train moves on a 
predetermined rail line with highly stable time 
schedule [T^Tq] from the origin station to its 
destination terminal. Thus the information of train 
trajectory and network resources can be obtained by 
the CC in advance with high accuracy (H. Liang and 
W. Zhuang, 2012; S. F. Xu etc., 2014). Each infestation 
with power P watts covers a segment of the rail line 
based on its wireless transmission range. Let T k and T k 
represent the time instants for the train to come into 
and go out of the coverage range of the hth infestation 
respectively, where h G J~C = {1,2, ... H}, and w.l.o.g., 
T{ = Tj and T^ +1 = T 0 . Thus, the trip duration is 
partitioned into H non-overlapped periods according 
to the entry time, i.e., [T^ , T^ +1 ] for all h. Then, each 
period is further partitioned into K h slots with fixed 
length T F , where K h = [(T k+1 — T k )/T F ] . Assume 
w.l.o.g. thatK h = K, Vh. When the train visits the hth 
infestation, the kth (k G K h = (1,2,... ,K}) slot begins 
and ends at times T k + (k - 1)T F and T k + kT F , 
respectively. 

The data of CM services is encoded and segmented 
into a large number of packets, and we denote the link 
capacity A hk as the maximum number of packets that 
can be delivered from the hth infestation to the VS at 
the kth slot, i.e.. 


\o, a n<t<v M , 


where t = T k + (k — 1)T F , Vk, h, and we assume that 
each packet has equal size of L p bits, and R(t) 
represents the maximum achievable data rate at slott. 
For a downlink OFDMA system with N subcarriers, 
the R(t) is showed by 


R(t) = £„ =1 Wlog 2 (l + 


Pn Ct)hg(t) \ 
N 0 W +ICI n (t)/ 


(2) 


where W denotes the bandwidth of the sub-carrier, N 0 
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represents the single-sided noise power spectral 
density (PSD) of each sub -carrier, p n (t) and h n (t) are 
the transmission power and channel gain of the nth 
subcarrier at thetth slot, respectively, and ICI n (t) is the 
inter-carrier-interference experienced on the nth 
sub carrier at thetth slot. 


The wireless channel is assumed to be frequency 
selective fading, where different sub-carriers will 
experience different channel gains, and the knowledge 
of channel state information is available at CC. Thus, 
the total power P can be allocated among the 
sub carriers by using water filling such that the rate can 
be maximized. Noticing that the ICI power ICI n (t) 
caused by Doppler shift is not coordinated among 
different sub carriers, but its average impact is 
considered. By doing so, the wireless system can adapt 
to the varying channel timely and easily. According to 
(Y. Li and L. Cimini, 2001; T. Wang etc., 2006), the 
ICI n (t) can be expressed as 


ICL (t) 


C fdTsY 


ViV 

*-n = l,l*n 


vi 0) 

(l—n) 2 } 


(3) 


where T s is the OFDM A symbol duration, and 
f d = vf c /c represents maximum Doppler shift with v 
the moving speed, f c the carrier frequency and c the 
velocity of light. A tight universal upper bound on the 
ICI power (Y. Li and L. Cimini, 2001) can be employed 
to ease the challenge in this paper 

ICIupperbound < ^ (2nf d T s yfP) 2 . (4) 

Consider a single trip of a train from an origin station 
to a destination terminal. An on-demand CM service 
s E S = {1,2, ...,S} is requested at time G s , and if 
Q s E Z + packets are delivered to the VS beforetheir 
playback deadline D s , a reward w s will be obtained. 
Moreover, for a slotted transmission, the small 
difference between G s (D s )and the beginning time of 
its generated slot is omitted. Since the data packets 
delivered before their playback deadline can be 
successfully decoded at passenger devices, the 
required minimum average service rate to play out a 
services smoothly is showed by V s = Q S /(D S — G s ). In 
the sequel, we will thus investigate the optimal 
resource allocation for CM service scheduling. 


Problem Formulation 

In this section, we discuss an optimal resource 
allocation problem in the HSR downlink OF DMA 
system. The problem is formulated as a two-stage 
optimization programming. 

Define x hks as the packet number of service s delivered 


from the hth infos tat ion at kth slot. Let ^ be a delivery 
indicator of service, which equals to 1 if service is 
completely delivered before its deadline, and 0 
otherwise. Then i[> s can be expressed by 

{ 1, if ^2 Xllks - ® * ■ 

h€Uk€K. h 

0, otherwise. ( 5 ) 

The first stage optimization problem with the goal of 
maximizing the total reward of delivered services over 
a trip of the train is then formulated as 

(PI) max w sV s (6a) 

{Xhks} ^ 

s.t. Xhks € Z + , Vk 6 JCh, h eH.s E S (66) 

^2 ^3 Xhks - € s ( 6c ) 

h£H k£lC h 

y j Xhks < Ahk , v/b E /Cfc, h EH (6d) 

ses 

%hks — 6, if G s > Tf, + kTy 

or D s < Th T (k — 1 )Tp, 

V/c E JCh , h EH, s E S (6e)(6) 

where the constraint (6b) means that negative resource 
allocation is not allowed, (6d) states that the packet 
number that can be delivered from hth infestation at 
k th slot is limited by the capacity A hk , while (6e) 
shows that the packets can only be delivered between 
the request time and the deadline. 

The problem PI will determine the optimal service set 
(with maximum total reward), denoted by S* , that can 
be completely delivered before their deadlines. Notice 
that S * might not be unique. And moreover, PI can 
only guarantee the data transmission integrity, 
without considering of the continuity of data streams 
as well as the arrival order. As we known that because 
data delivery over the wireless links with limited 
resource and intermittent connectivity, the strict 
continuity constraint is prohibitively difficult to meet. 
But on the other side, due to the human perceptual 
inability to notice small deviations of presentation 
rate, people are more sensitive to a low-frequency but 
long-lasting disruptions than to a high-frequency but 
short-lived disruptions (M. C. Yuang etc., 1996). As a 
result, a better perceptual quality can be expected by 
replacing large continuity disruptions with shorter 
ones (N. Laoutaris and I. Stavrakakis, 2002), and it is 
beneficial to minimize the amount of cumulative 
discontinuous packets during each slot to improve the 
perceptual quality. The cumulative discontinuous 
means that the amount of available data packets 
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within each slot is less than that consumed by 
playback at the same time. Accordingly, the 
cumulative discontinuous packet number of service at 
the k th slot when the train is served by the hth 
infestation can be formulated by 

( h-l K h-l k 

Z Z + Z 

m = 1 t = 1 t= 1 

(7) 

whereA t hks = [min { T k + (k - 1)T F , D s } - G s ] denotes 
the transmission duration of the service from G s to the 
kth slot within the hth infestation, V s A t hks and 
Z^"! 1 ! Zf^ 1 x mts + Z t k =1 x hts refer to the required packet 
number for continuous playback of multimedia data 
and actual transmitted packet number during A t hks for 
service, respectively. 

Then, the second stage problem with the objective of 
minimizing the weighted total number of cumulative 
discontinuity packets of the delivered services in S* 
over a trip of the train is formulated by 

( p2 ) min ^w„f(x hk „) (8a) 

hen keKh ses- 

s.t. x hks € Z + , Vk e K h , h € n, fi e S\ (8b) 

E E Xhk* < Qsi Vs eS*. (8c) 

he W keic.h 

y X hks < A hk , Vfc e !C h ,\/h e H, (8 d) 

* es* 

x hks = 0, if G s >Ti + kT F 

or D s < Tj t -\- (h — 1 ) Ty , 

Vk € 1C h , V/i. € H,Vs € 5* . (8c)( 8 ) 

The potential rational under the above formulation is 
to facilitate data transmitting, and if a data packet is 
late, it should be transmitted as early as possible, 
instead of being discarded, to protect the continuity of 
the stream from further degradation. Moreover, 
services with higher reward values receive higher 
priority to obtain enough packets for maintaining 
continuous playback. 

Offline Resource Allocation Algorithm 

As we known, two- stage optimization is, in general, 
significantly harder than single-stage optimization. In 
this section, we will reformulate the two-stage 
problem P1-P2 into an equivalent single-stage 
optimization problem by fully exploiting the inherent 
structure. 

2-in-l Reformulation 

To proceed, let us define the weighted reward matrix 


W as efwi, w 2 , • • • , wsl, where e G R HKh represents 
the all-one vector. The packet allocation matrix is 
denoted by X G R //k ^ XlS ’ with its entry[X] hkb = x h k h * 
. Then we have 

E E E (E E ^ + E *4 = |,r (w x BX). 

Where ^ is a well-defined (HK^ x II K k )~ dimensional 
unit fewer triangular matrix. Meanwhile, we note that 
Ylhen YlkeJCb S$eS w sVsAthks is a constant. Thus the 
problem P2 is equivalent to 

max tr(W T BX) 

s.t. (8b) - (8e). ( 9 ) 

Introduced a parameter <y for the two- stage 
programming, then we have 

(P3) max Ws'&s + atr(W T BX) 

s.t. (6b) - (6e), ( 10 ) 

Where < *■ is a constant satisfying b < a < ,. r ( W T B x* ) , 
w* = ■ -ws}r an ^ X" is the solution to 

problem P3. Actually, if there are multiple solutions to 
PI, the above formulation is capable of picking the one 
with minimum weighted total cumulative 

discontinuity packets. We can draw the following 
equality between P1-P2 and P3. 

Lemma 1. The two-stage optimization problem P1-P2 
is equivalent to the problem P3. 

Proof. The proof of Lemma 1 follows (Y. F. Liu etc., 
2012) and is omitted for brevity. 

The problem P3 is a mixed integer nonlinear 
programming (MINLP) problem, which in general 
NP-hard (L. A. Wolsey, 1998). Now, we can transform 
the single-stage problem P3 to an -norm 
minimization problem. Let cj) s = 1 — and we have 

0s = \\Qs - Z Zhen Sfce/Cfc x ^llo * 

Thus the problem P3 is equivalent to the following 
sparse optimization problem 


- atr(W T BX) 

, (11) 

s.t. (6b) - (6e). 

Now, it can be clearly seen that the difficulty of (P4) 
comes from the optimization over integer variables 
and the non-convex objective function involving t 0 - 
minimization, which is NP-hard in general (L. A. 
Wolsey, 1998). Therefore, a sub optimal algorithm is 


(P4) min ^w s 


Q* E E 


%h.k.8 


hen keJCi, 
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proposed in the following subsection to solve this 
problem efficiently. 


Linear Programming Approximate (LPA) 

First of all, since f 0 - optimization problem P4 is NP- 
hard, it is natural to consider its /^-convex relaxation. 
By relaxing the hard constraints x hks E lA to linear 
inequality constraints 0 < x hks < Q s , the final 
approximation of the problem P4 can be expressed as 


min w, Q., - yy x,a:J - atr(W T BX) 

{***.} sg5 h€Hk€Kk ! 

s.t. 0 < xt„ks < Q.s : V/r e JCh,\/h E Ltys E S , 

(6c) - (6c). 


(12) 


However, it can be easily verified that the case 
H hew^kex h x hks > Qs doesn't happen, because the 
resources are underutilized by delivering more than 
Q s packets for service 5. Therefore, (12) is equivalent to 
the following LP problem 

min y'wsf Qs - y" x hk8 ] -atr(W T BX) 

y henkeKh ) 

s.t. 0 <x hkx <Q,,\/kelCh,Vhe'H.yseS. ( 13 ) 
(6c) - (6e), 

which is a convex optimization problem and the 
optimum solution x kks can be obtained by CVX. Since 
x hks not necessarily a nonnegative integer, we need 
more steps to attain a feasible solution. First, the 
solution x kks is rounded down to the nearest integer 
x* ks = L x hks J * Then, based on the insight discussed in 
section 3, a packet loading scheme is detailed in 
Algorithm 1. 

In step 1, the remaining packet number of services 
and the remaining capacity of the kth slot from the hth 
infestation is showed respectively by 

Qs = Qs “ ^hejf^ke^h *hks ' (14) 

Ahk = A hk - ^ X* ks . (0.0.15) 

s ES 


In step 5, §hk represents the active service set which 
can possibly be allocated with more packets at the/cth 
slot of the hth infestation, showed by 

= | 5 | s £ $ Qs > Oi + kTp, D s >T 1 h -\-(k- 1)Tf| (16) 


Algorithm 1 iteratively allocates the remaining 
capacity A hk to the active services in descending order 
of w s /Q s , until the capacity of the slot is fully utilized. 
Based on these principles, the complete description of 
the proposed linear programming approximate (LPA) 
algorithm is showed in Algorithm 2. 


Algorithm 1 Packet Loading Scheme 
l: Initialize Q &7 x*^ a 7 Vh e tC ht se S ; 

2: for h = 1 to H do 
3: ibr h — 1 to Kh do 

4: while Ahk > 0 do 

5; Calculate by (16); 

6: if S* k = 0 then 

7: break; 

S: else 

9: s* = argmax{w s /Q^}; 

&€Sa 

11: := Qs* 

A-hk Ahk ~ L 

13; end if 

14: end while 

15: end for 

16: end for 


Algorithm 2 Proposed LP Approximate Algorithm 
1 : Initialize P,N,W,N 0 ,T S , fd,h n (t),Q s ,G s ,D s and w s , 
Vt E [Tj,To],Vft E H, k g K h ,s G 5; 

2: Calculate {Ahk} by water-filling; 

3: Resource allocation by solving (13). 

4: Let x* hks = an d call Alg. 1 for packet loading. 


Online Resource Allocation Algorithm 

In this section, we will study the online version of the 
offline algorithm that presented in the section 4. 
Supposing that the knowledge of the future service 
arrivals is not available, and only channel information 
per slot within a channel coherence interval is known 
at infestations. The main idea behind our online 
algorithm is to allocate the packets according to the 
offline policy based on the current service arrivals and 
channel capacities within the channel coherence 
interval, and then reallocate when new services arrive. 

Let L represent the number of slots in a channel 
coherence interval. Assume that a set S 0 packets arrive 
at time instant 1 , and meanwhile L slots channel 
capacities are available. Then we construct a (sub-) 
system over the slots [1, L]. For such a subsystem, we 
run the proposed offline algorithm to find the optimal 
packet schedule policy. The packet delivery follows 
this policy until new services arrive at t. Then we treat 
the current time instant t as new " 1 " instant, and 
update the set S 0 including the newly arriving 
services, and removing the old services whose 
deadlines are past or packets all delivered. For other 
old services, we update their D s and Q s as D s — t and 
the number of their remaining packets, respectively. 
Provided by the new L slots channel information, the 
offline algorithm is runnig for the new subsystem 
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\t,t + L — 1] , and the subsequent packet delivery 
follows the resultant new policy. This process 
continues until all the packets are delivered. 

Simulation Results 

In this section, we evaluate the performance of the 
proposed resource allocation algorithms. The 
cellular /infos tati on integrated network under 

simulation assumes that there are H = 5 infostations, 
each with the transmission range 100 meters. An HSR 
train moves at the speed of 360 km/h and all wireless 
links experience Rician fading, where the channel gain 
h n (t) is generated by Rician distribution with K- factor 
1. The duration of a time slot is 20 ms. The data packet 
size is L p = 500 bits. Moreover, the size of each service 
(Q s ) is uniformly distributed within [25, 50] kbits. The 
service requests follow the Poisson process with 
average service arrival rate A, and the lifetime ( D s — G s ) 
of each service is exponentially distributed with 
average value 10 seconds. The reward of each service 
(w s ) is uniformly distributed within [1,10 ]. The other 
simulation parameters are set as below: T s = 1.33 x 
10 -4 s, f c = 3.5 GHz, W = 7.5 KHz, N 0 = 2sl0~ 7 
W/Hz, P = 20 dBm, L = 4, N = 4, e = 0.1, and the 
maximum number of iterations n max = 10 . 

We compare the performance of the proposed offline 
and online algorithms with three existing algorithms: 
first-in-first-out (FIFO), earliest due date (EDD) and 
the Smith ratio based algorithm in (H. Liang and W. 
Zhuang, 2012). All figures are obtained by averaging 
over 500 independent simulation runs. 

Fig. 2 illustrates the relationship between the total 
reward of delivered services over the trip of the train 
and the average service arrival rate (A).It can be seen 
that, by incorporating the demands of multiple CM 
services and train trajectory, our proposed offline 
resource allocation algorithm achieves the highest 
total reward.lt is also interesting to note that the 
proposed online algorithm performs close to the 
offline algorithm and outperforms the other three 
algorithms. Moreover, as shown in Fig. 3, the weighted 
total number of cumulative discontinuity packets 
increases as A grows from 0.5 to 1.2 . For a larger 
number services, less resources can be allocated to 
deliver each service. As a result, a less number of 
services can be delivered smoothly under the limited 
resources. It can be observed that the cumulative 
discontinuity packets of the existing algorithms, which 
don't take the continuity of data streams into 
consideration, are much more than those of the 


proposed algorithms. That means the proposed 
algorithms achieve better perceptual quality than the 
existed algorithms. 




FIG. 3 WEIGHTED TOTAL NUMBER OF CUMULATIVE 
D B CONTINUITY PACKETS VERS IB A 


From Fig. 2 and 3, we can see that the proposed 
resource allocation algorithms can obtain more 
rewards but less cumulative discontinuity packets 
than the existing algorithms under different system 
parameters. 

Conclusions 

In this paper, we investigate the optimal resource 
allocation for on-demand continuous multimedia 
service delivery in a cellular/infostation integrated 
network for high-speed railways. A two-stage 
optimization problem is formulated by taking account 
of the integrity and continuity of the on-demand data 
delivery. To address the NP-hard mixed integer 
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programming, an LP approximation algorithm has 
been developed. Guided by the optimal offline 
strategy, we proposed an online algorithm for 
practical systems where a-priori knowledge of the 
future service arrivals is not available. Finally, the 
simulation results demonstrate the effectiveness of the 
proposed algorithms. 
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